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Abstract 

Networks of steps, seen in STM observations of vicinal surfaces on Au and Pt 
(110), are analyzed. A simple model is introduced for the calculation of the free 
energy of the networks as function of the slope parameters, valid at low step 
densities. It predicts that the networks are unstable, or at least metastable, 
against faceting and gives an equilibrium crystal shape with sharp edges either 
between the (110) facet and rounded regions or between two rounded regions. 
Experimental observations of the equilibrium shapes of Au or Pt crystals at 
sufficiently low temperatures, i.e. below the deconstruction temperature of the 
(110) facet, could check the validity of these predictions. 

INTRODUCTION 

(110) surfaces of fee metals have been intensively studied in the past years by 
means of several experimental techniques. It was found by scattering experiments, and 
subsequently confirmed by direct STM observations that some "heavy" metals such as 
Au, Pt, Ir form a (2 x 1) reconstructed state. Other metals, e.g. Pb and Al, are not 
reconstructed. 

The (2x1) structure is also known as missing-row reconstruction because one 
out of two rows of atoms, aligned along the [110] direction, is missing from the surface 
layer. We stress also that there are two different realizations of these reconstructed 
states, in which either the even or the odd rows are missing. 




Figure 1: (a) Profile and (b) top view of a missing row reconstructed (110) surface. 

Fig. |l] shows a top and a side view of the missing-row reconstructed structure. 
Along the [001] direction the surface assumes a hill- and- valley profile, where the sides 
of the hills are actually (111) microfacets. In the (111) orientation, surface atoms are 
closely packed, therefore such orientations are energetically favored. 

In several STM investigations (Gritsch et al. 1991; Gimzewski et al. 1992; 
Kuipers, 1994) of vicinal orientations of (110) Au and Pt surfaces an unusual pattern 
of steps, as shown in Fig. [2], was observed. The figure represents a surface orientation 
slightly tilted towards the [110] direction with respect to the (110) facet. The missing 
rows, which are not visible in the figure, run along the vertical direction. 




Figure 2: STM images of vicinals of Au(lW) (courtesy of M.S. Hoogeman, L. Kuipers and 
J.W.M. Frenken, AMOLF Amsterdam). The area shown is of 190nm x 120nm at T = 550K, 
with a miscut angle of 0.07 degrees. 

This pattern of steps is unusual indeed. In normal situations a miscut along the 
[110] direction is generated by steps that run perpendicular to the missing rows, i.e. 
along the horizontal direction of Fig. [3|, because such an arrangement minimizes the 
total length and hence the total free energy of the steps required to produce the miscut. 
Instead, on missing-row reconstructed surfaces the steps zig-zag and repeatedly touch 
each other at a collection of contact points; one can also say that they form a network of 
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Figure 3: Examples of clockwise (CW) and anticlockwise (ACW) steps. 



two arrays of roughly parallel steps crossing each other and forming on average angles 
(f) and — </>, with the vertical direction in Fig. [2|. An explanation for the formation 
of this network has been given already in (Kuipers, 1994). In this article we give a 
more quantitative description based on a simple, yet, we think, quite realistic model, 
which describes surface orientations close to the (110) facet. On the basis of this 
model we discuss the thermo dynamical properties of the networks (Carlon and Van 
Beijeren, 1996). We find that the network is actually unstable or at least metastable: 
in equilibrium it decays into a combination of stable surface orientations, with the 
appearance of sharp edges in between. As a consequence the edge of the (110) facet 
shows cusps, i.e. jumps in the direction of its tangent. 



THE MODEL 

In a missing-row reconstructed (110) facet one can distinguish two different types 
of steps parallel to the missing rows, commonly known as clockwise and anticlockwise 
steps and illustrated in Fig. [3|. In STM observations anticlokwise steps are rarely seen 
especially at low temperatures, and when present they are mostly pinned by impurities. 
This indicates that clockwise steps have markably lower free energies per unit of length 
than the anticlockwise ones; this energy difference can in principle be derived from 
STM observations, by estimating the relative frequencies at which step segments of 
given length of either type of step do occur. 

On the basis of STM observations on Au(llO) surfaces it has been concluded 
(Kuipers, 1994) that at room temperature and for sufficiently clean samples, anticlock- 
wise steps should be absent from the surface. 

Due to the presence of reconstruction not all possible configurations of clockwise 
steps will actually occur on the surface. As illustrated in Fig. ^|(a) a clockwise step going 
up followed by another clockwise step going down induces a shift in the reconstruction 
of the lower terrace. For this reason, a closed clockwise step cannot be formed on 
a given terrace unless it is accompanied by domain boundaries separating regions of 
different reconstruction order (Fig. |3](b)). 

This fact may have some important consequences, for instance, on how the (110) 
facet roughens. In the present work however we will focus on the behavior of the system 
far from the deconstruction or roughening critical points. 

The interplay between step orientation and surface reconstruction is essential for 
understanding the formation of the network of steps. Indeed, a step perpendicular to 
the missing rows cannot zig-zag forming clockwise segments parallel to the missing 
rows, if expensive domain boundaries between opposite reconstruction states are to be 




(b) 

Figure 4: (a) A clockwise step up followed by a clockwise step down generate a shift in 
the reconstruction order; i?_ and R + indicate the two possible reconstructed phases, (b) As 
a consequence of this a closed clockwise step must be accompanied by domain boundaries 
between the two phases R- and R + . 
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Figure 5: Example of microscopic configuration of a step; if anticlockwise steps are to be 
avoided only segments along the +y and ±x directions are allowed. 



avoided. 

The single step free energy 

To clarify and quantify the previous discussion we calculate the free energy of an 
isolated step as function of the step orientation. For convenience we consider, for the 
energy of an anticlockwise step segment, the limit E^cw —> 00 : on ly clockwise step 
segments are allowed and they have energies per unit lengths 5 y and S x where y and x 
are the directions parallel respectively perpendicular to the missing rows. One could 
eventually introduce more parameters, such as a corner energy, but this is not essential. 

Let us consider the step as an oriented walk from a lattice point (0, 0) to (L x , L y ), 
as shown in Fig. |5|. As seen above, if anticlockwise steps are to be avoided, only step 
segments in the ±x-directions and the +y-direction are allowed. A simple calculation 
(for more details see Appendix I) gives the following free energy of a step per unit of 
length, tilted over an angle <ft with respect to the missing rows: 
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with A and B some non-negative constants. The second term on the 
the energy of segments of steps parallel to the ^-direction, while the third term is trie 
decrease in free energy due to entropy. The entropic term dominates at sufficiently 
small |e|, except at T = 0, where B = 0. 

Obviously f s {4>)/\ sin0| diverges for 0^0. It is minimal for an angle 0o satisfy- 
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by applying a one-dimensional 
orientation (Van Beijeren and 
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1 The quantity under the square root in (||) becomes negative at small values of f35 x and @6 y ; this 
happens at temperatures above the roughening temperature of the (110) surface when the solid-on-solid 
approximation for the step free energy f s (4>) is not positive definite and the simple theory considered 
here breaks down. 
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Figure 8: Schematic view of the network of crossing steps. 

built out of an alternating sequence of segments of orientations <p and — O - Then at 
any non-zero temperature the lowering of free energy through entropy gain outwins 
the increase through the excess energy at the corners between the successive segments, 
provided the density of these corners is sufficiently low. This is a variant of the usual 
argument against phase transitions in one-dimensional systems at finite temperatures. 
In the present case this argument does not work since, due to topological constraints 
induced by the reconstruction, a configuration of short "zig-zagging" step segments of 
orientations O an d — <fio joined together cannot be formed, as pointed out above. 

The surface free energy 

Next we calculate the surface free energy as function of the slope parameters p 
and q, which denote the tangents of the tilting angles of the surface with the [110] and 
the [001] directions. Both parameters q and p will be considered small, so we restrict 
ourselves to vicinal surface orientations. 

Interactions between steps are introduced in a simple way: a free energy e is 
associated to each pair of crossing steps, e keeps account of both short and long range 
interactions between the steps and it is defined as the excess free energy associated to 
the intersection point of two crossing steps forming angles 4>i and 2 with the ^/-direction 
on a (110) facet: 

e(0!,0 2 ) = F^^-Afo-sJM-sMfo) (6) 

Here F(<f>i, fa) is the free energy of the surface of area A with two crossing steps of 
lengths Si and s 2 , f is the free energy of the (110) facet per unit of area (no steps 
are present), and f s {4>) the free energy of a single step, whose properties have been 
discussed above. 

One could try to estimate e from a microscopic point of view taking into account 
long range interactions of elastic or entropic type and short range interactions near 
the crossing point. This estimate is non-trivial and for the general features of the 
equlibrium crystal shape that we want to address here, only the sign of e matters. In 



general, it is expected (Bernasconi and Tosatti, 1993) that elastic interactions between 
two antiparallel steps (as are the crossing steps of the network, provided O is n °t too 
large) yield a negative contribution to e. Entropic interactions (Gruber and Mullins, 
1967; Jayaprakash et al., 1984), which become relevant at higher temperatures, are 
repulsive and give a positive contribution to e. Finally at the crossing point one would 
expect a positive contribution to e. In the rest of the article we discuss both possibilities 
for the sign of e, which yield two different scenarios for the equilibrium shape of the 
crystal. 

Let us consider a miscut along the [110] direction, generated by a pattern of cross- 
ing steps forming angles 0i and 2 with the y-axis, as shown in Fig. The dashed 
inclined lines indicate hypothetical parallel isolated steps that would generate the same 
macroscopic orientation as the network; 1/q and 1/p are the average distances between 
these steps along the x and the y direction. 

The free energy per unit of projected area can be written as: 

f(p, q, X , 2 ) = f + ± [hf.fa) + l 2 f s ((f> 2 ) + e(0 1; 2 )] (7) 

where fo is the free energy per unit area of the (110) facet, while l\ and l 2 are the 
lengths of the two sides of the dashed triangle, of area A, shown in Fig. [8|. 

The actual free energy can be found by minimizing the free energy ([?]) with respect 
to the angles 0i and 2 ; to lowest orders in p and q and for ip < it/2— o , this amounts to 
minimizing f s (<f>)/\ sin0|. Therefore the minimum of (0) is at <pi w — 2 = 0o + 0(p, q), 
where 0o is the angle given by (J5|). Substituting this back into ([/]), expressing A, li, l 2 
as function of p, q, 0i, 2 and expanding to lowest orders in p and q we find: 

f(p,q) = /o + ^P+|fr^--tan0 o g 2 ) (8) 
sin 0o 2 \ tan 0o / 

where e is the interaction free energy (|) calculated at angles 0i = — 2 = 0o- As 
in usual expansions of surface free energies around facets (see e.g. Van Beijeren and 
Nolden, 1987) the term linear in p represents the contribution of non-interacting steps. 
The interaction terms are quadratic in the step densities, differently from usual step- 
step interactions which lead to terms cubic in the step densities (Gruber and Mullins, 
1967). The origin of the quadratic term can be understood easily: the number of step 
crossings per unit area is simply proportional to the product of the densities, p±q tan O , 
of the two types of steps. In addition there are long range interactions between the 
parallel steps in the network, but they will only contribute to terms of cubic or higher 
order in the step density expansion of the free energy (Gruber and Mullins, 1967). 

For ip > 7r/2 — 0o the expression (|7|) is minimized by a single array of steps (so 
0i = and l 2 = 0) and the free energy takes the usual form (Gruber and Mullins, 
1967): 

fM = fo + fs(<P)^P 2 + q 2 + 0((p 2 + q 2 )l) (9) 

Notice that the expression @ for the free energy of the step network, irrespectively 
of the sign of e, is a non-convex function of p and q. This result implies that the network 
is always unstable: some surface orientations disappear from the equilibrium shape of 
the crystal and are replaced by sharp edges. 



EQUILIBRIUM CRYSTAL SHAPES 



Figure 9: (a) The network of crossing steps, and (b) a combination of two step arrays of 
different orientation, after the phase separation, predicted by (5) for e > 0, has occurred. The 
+ and — indicate the opposite reconstruction phases induced by the clockwise steps. 

As is well known (see, for instance, Van Beijeren and Nolden, 1987), the equilibrium 
crystal shape is the shape that minimizes the total surface free energy at a given fixed 
volume. From the minimization of the free energy calculated above we can construct the 
equilibrium shape of the crystal around the (110) facet. This shape depends crucially 
on the sign of the interaction free energy e, therefore in the rest of the paper we will 
distinguish two different cases. 

e > 

When e is positive the system can decrease the surface free energy reducing the 
density of crossings. By applying the Maxwell construction to Eqs. (||, ||) one finds that 
a given surface orientation (p, 0) separates into two orientations of slopes (p,ptan0 o ) 
and (p, — ptan0 o ) joining at a ridge, as depicted in Fig. §(b). In this combination of 
two surface orientations all the crossing points have been "eliminated" with a net gain 
of free energy. To generate the interface between the two surface orientations one has to 
pay a positive amount of boundary free energy, but this free energy will be proportional 
to the length of the boundary and therefore negligible in the thermodynamic limit, when 
compared to terms proportional to the surface area. 

The structure of the domain boundary is quite interesting as well. Two clockwise 
steps running under angles 7r/2 + (fi and n/2 — </> cannot be simply joined together, 
because this will cause a mismatch in the reconstruction order. Such joints must be 
accompanied by a domain boundary. 

As we saw already in the previous Section the model predicts that the network is 
unstable, suggesting that the pattern of steps observed in STM experiments should de- 
cay, after a sufficiently long time, into the combination of surface orientations depicted 
in Fig. |9](b) In the present case this time may be extremely long. In fact once the 
network has been formed, the process of disentangling it into stable orientations as the 
ones shown in Fig. §(b) , may require the investment of a large amount of free energy to 
go through very unfavorable states and it may be difficult to observe it experimentally 
without a careful long annealing of the surface. 

Thermodynamically the metastability encountered here is highly unusual. In e.g. 
a homogeneous gas-liquid system a free energy that is a concave function of density 
always leads to instability due to spinodal decomposition. In our system this is im- 
peded by topological constraints on the steps, requiring concerted mass transportation 
over relatively large distances for the decomposition of a network into stable surfaces. 
Therefore even a non-convex free energy can be metastable. For describing the surface 




Figure 10: Arrangement of the steps surrounding the (110) facet for e > in the stable 
phase (a), in the metastable phase (b) and for e < (c). 

free energy one can use Eq. (Q) again. A typical arrangement of steps around the (110) 
facet in this situation is shown in Fig. TO(b). In this case the shape profile along the y- 
axis of the crystal for vicinal orientations is expected to be of the type z(y) ~ (y — yo) 2 , 
due to the term proportional to p 2 in the surface free energy. A free energy expansion 
with a term cubic in the step density would produce a shape profile with an exponent 
3/2, the so-called Pokrovsky-Talapov exponent (see, for instance, Van Beijeren and 
Nolden, 1987). 

e < 

When a negative free energy is associated to a crossing point the network of steps 
tends to condense to maximize the density of crossings. On the other side, entropic 
repulsions (Gruber and Mullins, 1967) favor configurations where steps are far apart 
and tend to stabilize the network. 

The competitions between repulsions and attractions may give rise to the two 
following scenarios: 1) if attractions dominate for the whole range of surface orientations 
the network fully condenses and the stable orientations along the p direction are the 
(110) and (100) facets, which would be directly connected under a sharp edge; 2) if 
repulsions, possibly of entropic type become dominant, then the surface free energy 
turns concave at higher step densities and the (110) facet is connected under an angle 
with the rounded regions. 

As pointed out above, the surface free energy f(p,q) contains terms which are of 
cubic or higher order in the step density. A term cubic in p with a positive coefficient 
could well stabilize the free energy at not too small values of p. In both scenarios the 
(110) facet has sharp edges running roughly perpendicular to the missing-row direction. 

A simple model with a negative crossing energy (Carlon et al., 1996; Carlon, 1996) 
the staggered body-centered-solid-on-solid- model (BCSOS model), is discussed in some 
details in Appendix II. The model reproduces both scenarios 1) and 2), depending upon 
temperature. 

Finally in Fig. [7| the dashed lines show the truncation of the equilibrium shape of 
the (110) facet by ridges connecting the facet to rounded areas and Fig. |T0|(c) shows 
the expected arrangement of steps around the truncated facet. Notice that in this 
case there are sharp ridges between rounded regions covered by networks of steps and 
regions covered by non-crossing step arrays. 



CONCLUSION 



In the present article we introduced a simple model which describes networks of 
crossing steps, that have been observed in STM experiments on Au(llO) and Pt(llO) 
surfaces. The model shows that these networks are unstable (or at least metastable) 
with respect to faceting. 

From the calculation of the surface free energy for orientations close to the (110) 
facet the equilibrium shape for this facet has been derived. For repulsive interactions 
between two crossing steps (e > 0) the facet shape is rather unusual: it is cusped 
and elongated along the [110] direction. The cusps are connected to ridges separating 
rounded regions. If the interaction is attractive (e < 0) the (110) facet is still elongated 
with sharp edges roughly perpendicular to the missing-row direction. In this case there 
are much weaker cusps, between these edges and the smooth facet boundaries roughly 
parallel to the missing-row direction. 

The predicted shapes could be observed in experiments on equilibrium shapes of 
crystals with (110) missing-row reconstructed facets. 

Heyraud and Metois (1980) studied shapes of small gold crystals in thermal equilib- 
rium with their vapor; in their samples only the (111) and (100) facets were observed, 
since the range of temperatures investigated (T m 1000°C) is above the roughening 
temperature of the (110) facet. To observe some of the shapes described in this article, 
temperatures below the roughening and the deconstruction temperatures of the (110) 
facet should be considered. 

We believe that the problems connected to the metastability/instability of the net- 
works of steps deserve further experimental investigation, as well. If e is positive, as we 
think should be the case for gold and platinum crystals, one should be able, starting 
from the metastable network of steps, to observe a nucleation of arrays of parallel steps 
connected under a ridge. Probably this could be observed in practice, e.g. in STM 
experiments, only for sufficiently high temperatures and in very pure samples. 
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APPENDIX I: CALCULATION OF THE STEP FREE ENERGY 

For the calculation of the step free energy we consider first the partition functions 
G x and G y of straight segments of steps along the x and y direction. The segments 
have energies per unit length equal to 5 X and 5 y . In Fourier space one has: 

+0O 1 1 

GM = E e~ 2n ^ (e**- + e~^) = _ + x _ (10) 

GJk v ) = Y e -^Sy/V2 -ink y = _± (n) 

n=l ° x 

Notice that segments of steps along the negative y direction are not allowed. We recall 
thatthe factors \[2 and 2 are the minimal length of step segments along the x and y 
directions. The step is generated by all possible combinations of horizontal and vertical 



segments: 
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The partition function of a step of lengths L x and is obtained by Fourier transform: 
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with L = yLl + and <p = axctaxi(L x / L y ). The integral over k y can be rewritten as 
a contour integral in the variable z = e thy ; the integrand has a simple pole in: 



The partition function now becomes: 
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Pf(k x ,<P) 



sin 6 + 



'^-v^loga + G?^)) 



COS ( 



(14) 



(15) 



(16) 



The dots in ([15]) denote terms which are not relevant in the thermodynamic limit 
L — > oo. Using the saddle point approximation, one can evaluate (|I~5"D as: 
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where the step free energy is f s (<t>) = f(ik,cj)) and i/c the stationary point of f(k a 
Working this out one obtains Eq. [L| 



APPENDIX II: THE STAGGERED BCSOS MODEL 

Without going into details we present here some relevant results concerning the 
staggered BCSOS model, which, for certain values of the energy parameters describes 
networks of crossing steps with a negative crossing energy (for more details the reader 
may consult (Carlon et al., 1996; Carlon and Van Beijeren, 1996 II)). 

The model is applicable to surfaces of ionic crystals of bcc type, as for instance 
CsCl, and describes all surface orientations (tsl) with \t\ + \s\ < 1. It does not have the 
characteristic feature, responsible for the formation of networks in the present paper, 
that an up step followed by a down step changes a reconstruction order; instead the 
networks are induced by a large corner free energy, leading to steps with long segments 
in the principal lattice directions, combined with a negative value of e, which favors 
crossings of steps. However, the model does illustrate the two different scenario's for 
step condensation discussed above for e < 0. 

At T = the crystal has the shape of a truncated piramyd, with a top (100) facet 
and four side facets of (110) type, as shown in Fig. |TT|(a). A step on the (100) facet has 
an energy per unit of length e > and a kink energy ex 3> e: steps consist, especially 
at low temperatures, of elongated straight segments with few kinks. These segments 
are oriented parallel to the two axes 71 and 72 shown in Fig. |TT|(b): differently from 
(110) surfaces of fee metals these two directions are equivalent. Steps can cross each 




Figure 11: (a) At T = the equilibrium shape of the BCSOS model is a truncated pyramid, 
(b) Steps are composed of elongated segments with few corners since the energy of a corner 
is much larger then the energy of a straight segment of step ex ^> e. Crossings of steps, 
indicated by circles in the figure, are energetically favored, i.e. at each crossing there is a gain 
of energy of 2e. 

other, with a gain of energy of 2e at each crossing point. In this model, thus, the 
crossing energy is e = — 2e < 0. 

The surface free energy has been calculated as function of the slope parameters p, q 
in mean-field approximation, for all possible values of p and q, also beyond the vicinal 
orientations. To give an example, we restrict ourselves to considering only f(p), i.e. the 
surface free energy for surface orientations of the type (lpO), with < p < 1; these are 
all the surface orientations along the axis a of Fig. [TT|(a). The mean-field calculation 
gives a surface free energy that at low temperatures is not stable (Fig. |T2|(a)): step 
attractions, due to the energy gained at crossings, dominate in the whole range of 
surface orientations and the network condenses until the optimal density of steps is 
reached, which in the present case corresponds to the (110) facet. In the equilibrium 
shape the two facets (100) and (110) touch each other at a sharp edge. The free energy 
instability is of the same form as that predicted in Eq.(8), with a negative value of e. 

We stress that the free energy f(p) of Fig. |12|(a) becomes convex for sufficiently 
large values of p. This change of sign of the second derivative of f(p) corresponds to a 
transition from a reconstructed to a non-reconstructed state. In the present case this 
critical point is of limited relevance since it occurs in an unstable region, and there 
seem to be no topological restrictions impeding the decay of an unstable orientation 
into two facets. 

The same mean-field analysis at somewhat higher temperatures gives a surface free 
energy of the type shown in Fig. |I3|(b), which is also concave at low p, but convex at 
higher values of p. At low step densities, attractions still dominate, yielding a concave 
(i.e. unstable) free energy. Entropic repulsions (Gruber and Mullins, 1967) which 
become important at high temperatures and high step densities, stabilize the network. 
The Maxwell construction gives as coexisting surface orientations the (100) facet and 
a non-faceted orientation of slope p . These touch again at a sharp edge. 
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Figure 12: Surface free energy along the axis a of Fig. |Xl| . (a) At low temperatures the free 
energy is not stable for the whole range of orientations < p < 1; the dotted line denote the 
Maxwell construction, (b) At higher temperatures the range of orientations po ~ 0.4 < p < 1 
is now stable. 
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